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Abstract. The opacity- modified Eddington limit has 
been computed for hydrogen-deficient model atmospheres. 
The R Coronae Borealis (RCrB) stars are found to be lo- 
cated strikingly close to the limit, which suggests that the 
unknown trigger mechanism for their visual declines of 
the stars are instabilities in connection with the stars en- 
countering the Eddington limit in their evolution. It also 
points to a similarity between the eruptive behaviour of 
the RCrB stars and the Luminous Blue Variables (LBVs). 

Super-Eddington luminosities in hydrostatic model at- 
mospheres manifest themselves by the presence of gas 
pressure inversions. Such inversions are not an artifact of 
the assumption of hydrostatic equilibrium but can also be 
present in hydrodynamical model atmospheres. Only for 
very large mass loss rates hardly realized in supergiants 
will the inversions be removed. Instabilities may, however, 
still be present in such inversions, which is investigated for 
both H-rich and H-deficicnt late- type supergiant model at- 
mospheres. Dynamical instabilities may occur in surface 
ionization zones, which might lead to ejections of gas. A 
\ local, non-adiabatic, linear stability analysis reveals that 
■ sound waves can be amplified due to the strong radiative 
forces. However, despite the super-Eddington luminosities, 
the efficiency of the radiative instabilities is fairly low com- 
pared to for early-type stars with growth rates of 10~ 5 s _1 . 

Key words: Stars: model atmospheres - Instabilities - 
Stars: AGB and post-AGB - Stars: mass-loss - Stars: evo- 
lution - Stars: variables: R Coronae Borealis - Stars: vari- 
ables: Luminous Blue Variables 



1. The Eddington limit and the development of gas 
pressure inversions 

The Eddington limit corresponds to a situation where 
the radiative acceleration outwards equals the gravita- 



tional acceleration inwards. Eddington (1926) originally 
only considered electron scattering for the opacity (the 
classical Eddington limit) , but in an ionization zone of an 
important element the true opacity can be much higher. 
Therefore, the opacity-modified Eddington limit in stel- 
lar atmospheres can occur at significantly greater gravi- 
ties than the classical Eddington limit, as is evident for 
example from the calculations of Lamers & Fitzpatrick 
(1988), Gustafsson & Plez (1992) and Asplund & Gustafs- 
son (1996). 

If hydrostatic equilibrium is required, a stellar luminos- 
ity which locally exceeds the Eddington luminosity auto- 
matically forces the development of a gas pressure (P gas ) 
inversion, as seen from the equation of hydrostatic equi- 
librium: 



dftot dP gas dP ra( j 



= -9P- 



dr dr dr 

From rearranging the above expression one obtains 



(1) 



ldP e; 



= -See = -g(l - 5rad/.g) = -g [1 - r(r)] , (2) 



p dr 

with the radiative acceleration defined by: 
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(e.g Mihalas 1978). Here F v denotes the physical flux and 
k u the total mass extinction coefficient (with dimension 
cm 2 g _1 ). Obviously a positive P gas -gradient must occur 
when the Eddington limit is encountered if hydrostatic 
equilibrium is valid. 

It should be emphasized that a density inversion does 
not necessarily imply that the Eddington limit is exceeded. 
A density inversion predominantly occurs due to a chang- 
ing molecular weight in the ionization zone of a dominant 
species, while a P gas -inversion reflects a super-Eddington 
luminosity due to high opacity (Asplund et al. 1997a). 
Since both are related to ionization they may often oc- 
cur in the same atmospheric layers, which has probably 
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caused the confusion found in the literature (e.g. Maeder 
1989). 

It should be noted that P gas -inversions only seem pos- 
sible to develop in optically thick conditions. If the radia- 
tive acceleration exceeds gravity when lines dominate the 
opacity, the radiative force is highly unstable to perturba- 
tions induced by velocity gradients, since the spectral lines 
may then absorb unattenuated continuum flux. There- 
fore, rather than the development of a P gas -inversion, a 
stellar wind will be initiated, which is the case for radia- 
tively driven winds of hot stars. The aim of the present 
study is to investigate whether P gas -inversions in late- type 
stars, where the opacity instead is dominated by contin- 
uous opacity, are subject to a similar instability which 
would prevent their existence. Before considering possi- 
ble instabilities, in Sect. || the derived Eddington limit 
for H-deficient stars is presented and a possible connec- 
tion between the limit and the declines of the R Coronae 
Borealis (RCrB) stars is discussed. 

2. The Eddington limit and the RCrB stars 

2.1. The Eddington limit in the T e g-logg diagram 

The location of the Eddington limit in a T c s-logg dia- 
gram can be estimated by model atmosphere calculations 
(Lamers & Fitzpatrick 1988; Gustafsson & Plez 1992). 
It should be emphasized that these models consistently 
include the radiative acceleration in the equation of hy- 
drostatic equilibrium. Therefore, the location of the limit 
is not obtained by extrapolation towards lower gravities 
nor docs it imply zero density, as has sometimes been 
claimed in the literature (e.g. Humphreys & Davidson 
1994; Nieuwenhuijzen & de Jager 1995). An interesting 
conclusion from the above studies is that the theoretical 
opacity-modified Eddington limit coincides with the upper 
luminosity limit for stars (Lamers & Fitzpatrick 1988). It 
is therefore possible that the Eddington limit, possibly in 
connection with rotation (Langer 1997a, b; Owocki & Gay- 
ley 1997), prevents the evolution of stars into the super- 
Eddington regime by the existence of radiative instabili- 
ties. These may give rise to the eruptions of the Luminous 
Blue Variables (LBVs) (cf. Humphreys & Davidson 1994). 

Here a similar investigation but for the H-dcficicnt 
RCrB and related stars will be carried out; a preliminary 
report can be found in Asplund & Gustafsson (1996). The 
presence of P gas -inversions in such model atmospheres was 
noted by Asplund et al. (1997a), who speculated that it 
may be the trigger mechanism for the enigmatic visual 
declines of the RCrB stars. The variability of these stars 
is likely due to the formation of an obscuring dust cloud 
close to the stellar surface in the line-of-sight, but how 
the dust condensation can proceed despite the relatively 
high equilibrium temperatures has long been an enigma 
(cf. Clayton 1996). The presence of shocks seems to be a 
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Fig. 1. r as a function of optical depth for a sequence of mod- 
els with RCrB abundances at constant luminosity and mass. 
All models have C/He = 1%. The T eS = 8 000K model is 
only partly affected by convection in the inner layers, while 
the hottest model is not at all. Increased molecular absorp- 
tion causes the slight increase in V towards the surface for the 
coolest model 

very promising mechanism for the gas to cool sufficiently 
and thereby trigger dust formation (Woitke et al. 1996). 

The models used here are recently constructed line- 
blanketed model photospheres with typical abundances of 
R CrB stars (Asplund et al. 1997a). In Fig. [l] the variation 
of r with optical depth is displayed for a sequence of mod- 
els corresponding to constant luminosity and mass (e.g. 
L* = 8- 1O 3 L and M* = 0.8 M ). The super-Eddington 
luminosities only occur in the deep atmospheric layers 
where He I ionizes and, as expected, where P gas -inversions 
occur. Since the RCrB stars presumably evolve towards 
higher T e g (Kilkenny 1982), the stars will move from the 
sub-Eddington regime into having super-Eddington lumi- 
nosities, as is seen in Fig. || where the theoretical Ed- 
dington limit in the T e ff-logg diagram is shown. Com- 
pared with solar abundances the limit is shifted towards 
higher T c g due to the higher ionization potential of He 
compared with H, but also extends to slightly larger logg. 
The location of the Eddington limit for higher T c g has not 
yet been explored, since a self-consistent hydrodynamical 
treatment is needed for T c s 13 000 K when T > 1 occurs 
for tross < 1 and line opacity may determine g r ad- 

The exact location of the computed Eddington limit 
will naturally depend on the detailed properties of the 
model atmospheres. For T e g <, 8 500K convection in the 
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Fig. 2. The location of the opacity-modified Eddington limit (solid) in relation to the observed values of RCrB stars (Lambert 
et al. 1997; Asplund et al. 1997c) . Also shown is the classical Eddington limit (dashed) and the hot R CrB and EHe stars with 
determined parameters (Jeffery 1996). The cross in the upper left corner illustrates the typical uncertainties in the estimated 
parameters of the R CrB stars 



deeper layers of the He I ionization zone restricts the limit 
to lower logg, since the radiative flux diminishes in the 
presence of a significant convective flux. Also, the lower 
temperatures decrease k v and g ra d further. Including tur- 
bulent pressure would therefore bring the limit towards 
greater gravities by typically 0.2 dex, due to less efficient 
convection with lower densities. A C/He ratio of 1% has 
been assumed on the basis of the measured ratio in ex- 
treme helium (EHe) stars (Jeffery 1996), which presum- 
ably the RCrB stars are related to (Lambert et al. 1997). 
A higher carbon abundance pushes the Eddington limit to- 
wards higher gravities the same way as turbulent pressure 
does. A lower overall metallicity by 0.5 dex would decrease 
the limit by about 0.2 dex by lowering the temperatures 
in the inner layers with a diminished line-blanketing. 

2.2. Possible consequences for H-deficient stars 

In Fig. |2| the estimated parameters of the RCrB stars 
(Lambert et al. 1997; Asplund et al. 1997c) are also shown. 
Obviously, rather than being randomly located in the dia- 
gram as one a priori might have guessed, the R CrB stars 
fall close to the limit, which suggests that a connection 
between the declines of RCrB stars and the Eddington 
limit exists. Thus, there is an alluring resemblance be- 
tween the RCrB stars and the Luminous Blue Variables 
in their proximity to the Eddington limits and their erup- 
tive behaviour (Asplund & Gustafsson 1996). If the RCrB 



stars are evolving towards higher T c g (Kilkenny 1982) and 
the EHe domain at constant luminosity, their immediate 
progenitors should be the H-deficient carbon (HdC) stars 
(Lambert et al. 1997). The EHe and HdC stars experi- 
ence no visual fadings despite being otherwise similar to 
the RCrB stars. It is therefore plausible that the reason 
for the non-variability is that the HdC stars have not yet 
reached the Eddington limit while the EHe stars are lo- 
cated on the stable side of the limit at higher T e s. 

If the declines of the R CrB stars are indeed triggered 
by the super-Eddington luminosities, one could expect a 
correspondence between the maximum value of T in the 
atmosphere and frequency of declines. In the pulsation- 
induced dust condensation model for the RCrB stars 
(Woitke et al. 1996), in which the dust formation is trig- 
gered by propagating shock waves due to ordinary pulsa- 
tions, one would sooner expect a correlation with pulsa- 
tion strength. Intriguingly, V854Cen, which is the RCrB 
star most frequent in decline, has very weak photospheric 
pulsations (Lawson & Cottrell 1997) but by a wide mar- 
gin the largest T found of the 18 RCrB stars analysed 
sofar. RY Sgr, on the other hand, has by far the strongest 
photospheric pulsations but only average amount of de- 
clines, which might reflect the fact that T is not unusually 
high. Pulsation may, however, push some stars across the 
limit by periodically shifting them towards higher T e ff, e.g. 
AT off = ±250 K for RCrB itself (Rao & Lambert 1997). 



4 



M. Asplund: The stability of stellar atmospheres at the Eddington limit 



The proposed connection between the Eddington limit 
and the declines of the R CrB stars is certainly not without 
problems, which deserve to be addressed further. A few 
stars seem not to behave as expected if the proposal is 
correct, in particular this is true for some of the coolest 
EHe stars and the three hot R CrB stars. 



In the presence of a stellar outflow, the equation of 
hydrostatic equilibrium must be replaced with the mo- 
mentum equation: 



dv 
'dr" 



ldP„ 



dr 



ffoff, 



(5) 



3. The disappearance of gas pressure inversions 

3.1. Effects of convection 

As seen above, a P gas -inversion must develop at the Ed- 
dington limit in hydrostatic equilibrium. One should ask, 
however, whether strong convection may take place to pre- 
vent it, since the ionization zone causing the high opac- 
ity and thereby the super-Eddington luminosities also 
favours convection. In fact, a super-Eddington luminosity 
in the diffusion regime automatically requires the presence 
of convection as shown by Langer (1997a). Generalizing 
Langer's derivation for a non-negligible radiation pressure 
to also include the effects of ionization, the Schwarzschild 
criterion V ra d > V a( j for the onset of convection is equiv- 
alent to 

4a{2 + 8a + a;(l - x) [5/2 + e/fcP + 4a]} 

T(r) > LL - L '- — n (4) 

5 + 40a + 32a 2 + x (1 - x) [5/2 + e/kT + 4a] 2 

(e.g. Mihalas & Weibel Mihalas 1984). Here 
a = P ra d /P gas , while x denotes the ionization fraction of 
the most abundant element with an ionization potential 
e. Convection must be initiated as T(r) approaches unity, 
since the right hand side is always less than 1. Thereby 
the radiative flux and hence T diminishes when convec- 
tion is efficient. The Eddington limit can therefore not 
be exceeded in adiabatic convection zones in the stellar 
interior and hence P gas -inversions will not develop. Thus, 
inversions can only be present in non-adiabatic convection 
zones close to the stellar surface. 

It should be remembered that the location of the Ed- 
dington limit is dependent on the convection treatment. 
A Pg as -inversion will for example disappear if the den- 
sity scale height is applied instead of the normal pres- 
sure scale height for the calculation of the convective flux 
through the mixing length theory. The mixing length the- 
ory is rarely a good description of stellar convection, and 
in particular for these extreme situations it may be very 
misleading. Ultimately, hydrodynamical simulations will 
be necessary to investigate the effects of convection on 
Pgas-inversions. 

3. 2. Effects of mass loss 

Using the momentum equation one can investigate under 
which conditions a P gas -inversion can exist in the presence 
of a stellar mass loss. Or equivalently: which mass loss 
rates can one expect from a star if the super-Eddington 
luminosities would instead drive a stellar wind? 



where we have restricted ourselves to the one-dimensional, 
time-independent case. Also other external forces, such as 
turbulent forces (gtm-b = — l//9-dP tur b/dr) and centrifugal 
forces (<7ccnt = v rot( r i@)/ r )j can be contained in g e s be- 
sides radiation. It should therefore be kept in mind that 
the possibility of an instability may be underestimated 
when restricting the following discussion to only radiative 
forces, as other forces also tend to be destabilizing (e.g. 
Nieuwenhuijzen & de Jager 1995; Langer 1997a, b; Owocki 
k Gayley 1997). 

Together with a constant, spherically symmetric mass 
loss rate M = 4irr 2 pv, and the equation of state 
Pgas = TZpT/fJ,, it is possible to rewrite the momentum 
equation into the so called "wind equation": 



4) 



1 dv 

v dr 



1 dT 
T dr" 



1 dp 
p dr 



ffeff, 



(6) 



where the isothermal sound speed defined by c T = 
(<9P g as/<9p)T = Pgas/p has been introduced. By an in- 
spection of Eq. ||, it is possible to determine under which 
circumstances p- or P gas -inversions are allowed. In nor- 
mal cases, d/z/dr> in the surface layers and therefore 
the parenthesis on the right hand side of Eq. |^ is pos- 
itive as long as a temperature inversion is not present, 
which is not possible in optically thick layers with a ra- 
diative flux. Therefore, if g ra d > 9 but v < ct, the ve- 
locity gradient dw/dr must be negative. Since a constant 
M has been assumed, a density inversion will be present if 
r/v-dv/dr+2 < and therefore a P gas -inversion is possible 
(though not required). If, on the other hand, v > ct while 
flrad > 5; then dv/dr must be positive, and the density gra- 
dient cannot be positive: a density inversion cannot exist. 
Since a P gas -inversion requires a density inversion, neither 
inversion can occur in the stellar atmosphere despite the 
super-Eddington luminosities. Again the possibility of a 
temperature inversion has been neglected. A sufficient cri- 
terion for the disappearance of a P gas -inversion is therefore 
that the outflow velocity due to mass loss is high enough 
(v > ct). 

The critical mass loss rate can be estimated by com- 
paring v — M/4irr 2 p with or- The values for ct and p 
will be taken from static models, which is justified for 
these rough estimates, since the atmosphere below the 
sonic point is normally little affected by mass loss. A 
model atmosphere with T c g = 7 000K, logg = 0.5 [cgs] 
and solar abundances requires a critical mass loss rate 
M crit « 6- 1O -4 (AP/M ) Mq yr" 1 . Hence P gas -inversions 
can exist in stars despite high mass loss rates. Achmad & 
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Lamers (1997) have arrived at similar conclusions and es- 
timates. They also verify the finding numerically with dy- 
namical steady state atmospheres, though not altogether 
self-consistently. 

An even higher critical mass loss rate is found for 
the RCrB stars. The lack of hydrogen makes the con- 
tinuous opacity significantly lower than for solar abun- 
dances and hence the density at a given optical depth 
will be correspondingly greater, typically by a factor of 
50 (Asplund et al. 1997a). Since the density is higher, 
the outflow velocity will be lower in order to carry the 
same M. A similar model as above but with abundances 
typical for RCrB stars, will therefore require M cr it ~ 
3 • Kr 2 (A/*/M ) M yr" 1 . The mass loss rates of RCrB 
stars are poorly known, in particular since it may be 
episodic, but rough estimates suggest M ~ 10 -6 Mq yr -1 
(e.g. Feast 1986), orders of magnitude smaller than M cr it- 

It is concluded that P gas -inversions are not an artifact 
of hydrostatic equilibrium but may also exist in dynamical 
steady state atmospheres, even with a high mass loss rate. 

4. Dynamical instabilities 

Dynamical stability has been extensively studied, in par- 
ticular in connection with the termination of the AGB 
(e.g. Paczynski & Ziolkowski 1968; Wagenhuber & Weiss 
1994). Radial adiabatic oscillations will grow exponen- 
tially in time for dynamical instability, which is formally 
encountered if the first generalized adiabatic exponent 



ri = 



\ dhip 



(7) 



ad 



is less than 4/3 in a region effectively isolated from the 
rest of the star. An ideal gas has Tx = 5/3 while pure 
radiation has I\ = 4/3, but the combined effect of non- 
negligible radiation pressure and ionization may push Ti 
below the limiting value. For dynamical instability to oc- 
cur, the pressure weighted volumetric mean value of Ti 
between the depth r and the stellar surface P* 



<ri> = 



riPtotd^ , 



Ptotdy 



(8) 



must be reduced below 4/3 in a significant fraction of 
the star. Stars with a pronounced core-envelope structure 
and high luminosity-to-mass ratios, like the supergiants 
investigated here, may show dynamical instabilities in ion- 
ization zones (e.g. Stothers & Chin 1993; Wagenhuber & 
Weiss 1994). 

In Fig. U the run of Ti and (Ti) as functions of optical 
depth are shown for two model atmospheres correspond- 
ing to T e g = 7 000 and log g = 0.50 for a solar and a H- 
deficient composition. In particular, (ri) is reduced much 
below 4/3 for the model with solar abundances in the Hi 
and He I ionization zones. At the surface P ra( j dominates 
the total pressure and hence Ti is very close to 4/3. A 
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Fig. 3. a The generalized first adiabatic exponent Ti as a func- 
tion of optical depth in model atmospheres with solar (solid) 
and RCrB (dashed) compositions, b The pressure weighted 
first adiabatic exponent (ri) for the same models as in a. In 
both panels the parameters are T e g = 7 000 and logg = 0.50. 
Values below the dotted lines at Ti = (ri) = 4/3 correspond 
formally to dynamical instability 



more complete study for H-rich model atmospheres has 
been presented by Lobel et al. (1992), who arrive at the 
same conclusions. In the R CrB model only the He I ioniza- 
tion zone exists, which explains the quite different depth 
variation of Ti. The higher temperatures at depth in the 
H-deficient model make He I ionization occur at smaller 
tross • Also seen is the minor effect of the C I ionization 
zone at tr, OS s — 1- These models only extend down to 
tross = 40 and hence it is not possibly to tell how (Ti) 
varies further in, though the He II ionization zone will also 
reduce Ti. 

It seems that atmospheres of late-type supergiants are 
close to being dynamically unstable, or may even be so 
for sufficiently high luminosity-to-mass ratios. Violent in- 
stabilities due to this might lead to ejection of material, 
as found for the termination of the AGB (Wagenhuber & 
Weiss 1994), as well as for the LBVs (Stothers & Chin 
1993). Also RCrB stellar models seem to suffer from a 
similar dynamical instability as the stars on the tip of 
the AGB (A. Weiss, private communication). It should be 
remembered though that the above-mentioned evolution- 
ary studies lack a proper hydrodynamical treatment, and 
therefore interpreting the instabilities found in the models 
as actually occurring in the stars must still be made with 
some caution. 
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5. Radiative instabilities 

5.1. Radiatively driven sound waves in P gas -inversions 

It is clear from the discussion in Sect. [| and || that P gas - 
inversions may exist in stars with high luminosity-to-mass 
ratios and that such stars are close to being dynamically 
unstable. Thus, these stars may be vulnerable to radia- 
tive instabilities, since the additional destabilization from 
radiative forces may render them unstable. 

The existence of radiatively driven sound waves in 
hot stars due to high g ra d has been suggested by Hearn 
(1972, 1973), and further elaborated by Berthomieu et al. 
(1976), Spiegel (1976) and Carlberg (1980). The instabil- 
ity is based on the increase in opacity upon compression, 
since the radiation force is in phase with the velocity per- 
turbation and hence causes amplification. For these high 
temperatures it requires essentially isothermal perturba- 
tions, since 



(dK F /dp) T + (r 3 - l)(dK F /dT) p < 0. 



(9) 



Here T3 — 1 = (91nT/<91np) ac i, while kf is the flux weighted 
mean opacity. Thus only short wavenumbers will be am- 
plified. The growth rate is especially high as T approaches 
unity (Berthomieu et al. 1976) but the super-Eddington 
case has not been investigated. Carlberg (1980) found that 
the instability may provide rapid transfer of momentum 
from the radiation to the gas with typical growth rates of 
10 _3 s _1 for T cff ~ 50 000K. It is possible that a similar 
instability is exhibited also in late- type supergiants, with 
the important differences being the super-Eddington lumi- 
nosities in the P gas -inversions and that isothermal pertur- 
bations are not necessary, since both (<9kf/c?p)t > and 
(<9kf /dT) p > in the upper layers of the P gas -inversions 
(Fig. |). 

To study the stability of a P gas -inversion a local, lin- 
ear, non-adiabatic stability analysis has been carried out. 
Small, sinusoidal perturbations around the equilibrium 
values are assumed for the temperature, density and ver- 
tical velocity: 



i(kr-ujt) 



P = Po + Pi ■ e 
v = v Q +vi .e t(kr ~ ut \ 

where Ti, pi and v\ are small. The atmosphere is assumed 
to be static (vq = 0), since the outflow velocity must any- 
way be less than the isothermal sound speed for a Pgas- 
inversion to exist, as demonstrated above. The analysis 
will be restricted to vertical motion and therefore multi- 
dimensional instabilities of Rayleigh- Taylor-type will not 
be found. Carlberg (1980) investigated also such a gradient 
instability for hot stars but found slow growth rates com- 
pared to the radiatively driven sound waves. Such a gra- 
dient instability would resemble convective motion which, 
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Fig. 4. The dependence of the Rosseland mean opacity on 
temperature and density for a solar abundances and b RCrB 
abundances. The atomic data for the opacities are taken from 
the Opacity Project (Seaton et al. 1994 and references therein). 
The top of the P gas -inversion for a model atmosphere with 
T e ff = 7000K and logg = 0.5 occurs at logT « 3.9 for solar 
abundances and at logT « 4.1 for an RCrB composition 



however, is already present in the layers of interest in the 
unperturbed model and is thus not specifically searched 
for here. 

The equations to be linearized are the material equa- 
tions for conservation of mass, momentum and energy: 



P 



Dp 

Dt 
Dv 

Dt 
De 
Dt 



• P (V • v) 
1 



= 



VP gas = - 



Scffr 



Pgas Dp 



P 



Dt 



4 ^Ro S ss (J-B)-pv 3rad f 



(10) 

(11) 
(12) 



(e.g. Mihalas k Weibel Mihalas 1984), where D/Dt, e, B 
and J denote the Lagrangean co-moving derivative, the 
specific internal energy of the gas, the wavelength inte- 
grated Planck function (B = B u (T)dv) and mean in- 
tensity (J = Jq 00 J v dv), respectively. The absorption and 
Planck mean opacities have been replaced by the Rosse- 
land mean opacity for absorption k 



a bs 
Ross 



K Ro 



-a, which 



is a good approximation for the relevant optical depths. 
The scattering contribution to the total extinction is here 
denoted by a. Furthermore LTE (S„ = B v ) has been as- 
sumed. 
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The material equations must be coupled to the equa- 
tions for the radiation field to account for radiative trans- 
fer effects through exchanges of energy and momentum 
between the photons and the gas. In this exploratory in- 
vestigation the effects of radiation are included with a 
time-independent treatment, i.e. the radiation field is as- 
sumed to be quasi-static. The radiation momentum and 
energy equations will then look like: 



dK 

~ch 
dH 

~ch 



— pKpH 

P (KRoss -<J)(B - J). 



(13) 
(14) 



Here H and K denote the wavelength integrated Edding- 
ton flux (H = (l/47r) J °° F v dv) and radiation pressure 
(K = (c/47r)P ra d) respectively. In the diffusion regime for 
the unperturbed atmosphere, kf = kross is valid. The 
equilibrium radiation quantities are related by: 



J = B = 4(T /T eS ) 4 H = (T R I$/n. 



(15) 



For the first relation radiative equilibrium is implicitly im- 
posed, which is not strictly true in the cases which will be 
investigated below. However, the analysis will be restricted 
only to atmospheric depths where the convective flux car- 
ries < 10% of the total flux (i.e. in the upper part of the 
convection zone) and radiative equilibrium can therefore 
be largely justified. To close the relations the Eddington 
approximation (K = J/3) will be made, which is, like 
LTE, well justified in the atmospheric layers of interest. In 
the following, the radiation energy and momentum equa- 
tions will be combined to a single equation, which means 
that there is one more variable to perturb, for which J is 
chosen: 



J 



Jo + Ji -e 4(fer - wt) . 



Together with the definitions of the isothermal sound 
speed Ct, the density scale height H p = — [dlnp/dr]" 1 and 
an ideal gas for P gas , linearization to first order yields: 



^[-iu]+ Vl [ik-HZ 1 ] = (16) 
Po 



Pi 
Po 



Po 



To 



9 + Ct (ik + ^r^- - — ^ ] 
V T dr po dr J 

2 {■, „-i, 1 dT 2 dp\ 

r ■ l Jl 

+Vi [-IUJ\ + — 

•Jo 



+ 



+ 



-ik 



4crJq 



3cp 



+ 

= (17) 



Table 1. Local values of the atmospheric structure. Both the 
model with a solar and an R CrB composition have T c g — 7000 
and log g = 0.5 



Variable 


Units 


Solar 


RCrB 


log TRoss 




0.5 


1.2 


5rad 


[cms -2 ] 


1.29 • 10 1 


4.31 • 10° 


T 


[K] 


9610 


14950 


P 


[gem 3 ] 


4.62 • 10- 11 


4.96 • 10- 9 


p 


[dyn cm -2 ] 


5.19 • 10 1 


1.92 • 10 3 


Cs 


[cms -1 ] 


1.32 • 10 6 


6.67 ■ 10 5 


(l/T)dT/dr 


[cm -1 ] 


-6.67- 10- 12 


-3.15 ■ 10- 11 


(l/p)dp/dr 


[cm- 1 ] 


2.10 • 10- 11 


9.87 ■ 10- 11 


(l/Ai)<Wdr 


[cm- 1 ] 


5.92 • 10- 12 


6.50 • 10- 11 




[cmV 1 ] 


2.50 • 10° 


9.11 ■ 10- 1 


a 


Kg" 1 ] 


2.55 • 10- 1 


2.94 • 10- 2 


dlnn/dlnT 




2.43 • 10- 1 


1.32 ■ 10 1 


dlnn/dlnp 




6.75 • 10- 1 


1.72 ■ 10- 1 


de/dr 


[cm 1 s- 2 ] 


-1.29 • 10 2 


-5.15 • 10 2 


p(de/dp) T 


[cm 2 s- 2 ] 


-2.15 • 10 12 


-8.21 ■ 10 11 


T{de/dT) p 


[cm 2 s- 2 ] 


2.14 • 10 13 


1.63 • 10 13 



'To 



Pi 
Po 
de 



de 2 



iUjT °W + 16 ^° ~ Cr °' ) CrRT ° 4 

5rad + ^ + C T H p 



+Vl 

, Jl 



+ 
+ 



pi 

Po 



ik 



Jo 
1 fT cS 



[-4(K -<To)a R T 4 ] = (18) 



Ti 
'To 



Po(n - cr ) 4 V T 
ik 



^ 91n«\ 
dlnp J 



1 T 



dhiK 

po(no — cro) 4 V To ) d\vH 
k 2 



i + 



3k (ko - oo)pI 



= 0. (19) 



The notation has been abbreviated according to 
Pp = (Po/po)(dp/dp) and p T = (T / p )(dp/dT). Deriva- 
tives with respect to v and J have been ignored, which 
is justified in optically thick layers where line opacity is 
unimportant. True enough, photoionization edges intro- 
duce some velocity dependence on g rac j but this minor ef- 
fect is neglected. 

In order for the above set of equations to have a non- 
trivial solution for the perturbations Ti,p\, v\ and J 1; 
the determinant of the equations must be equal to 0. The 
three complex eigenvalues w for a given real wavenum- 
ber k then correspond to the allowed perturbations. Two 
of the modes will represent radiation-modified gravity- 
acoustic waves while the third is a non-propagating, heav- 
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ily damped, thermal wave. Here two different cases will 
be investigated numerically, one corresponding to a P gas - 
inversion in a supergiant with solar abundances and the 
second the same but for an RCrB star. As before, both 
model atmospheres have T c g = 7 000K and \ogg = 0.5. 
Values for some of the physical parameters at some atmo- 
spheric layer where the instability may occur are found in r ^ 1 
Table |. '„ 

The calculated growth rates are shown in Fig. pL which 
reveals amplification for k < 2 • 10 _8 cm _1 (solar com- 
position) and k < 2 • 10 _7 cm _1 (RCrB composition). 3 
For large k, the perturbations are as expected damped 
due to the negligible optical depths of the disturbances 
(ta = 2Trnp/k), which make the sound waves little af- 
fected by the radiative forces. Instead, the smoothing ef- 
fects from radiation energy exchange due to large gradi- 
ents in the perturbed variables dominate. The maximum 
growth rates are slightly smaller for a solar composition 
than for a H-deficient composition: 1.3 ■ respectively 
4.3 • 10~ 5 s -1 , which corresponds to an e- folding time on 
the order of the sound crossing time of the atmosphere. 
For the largest wavenumbers with such growth rates, the 
real parts are significantly higher, which implies that the 
sound waves are only slightly amplified per wavelength. 
For smaller k, however, the growth rate is comparable 
with the oscillation rate. In fact, for the model with so- 
lar composition, the wave is essentially non-propagating 
for the smallest k. However, the use of a local analysis is 
not justified when the wavelengths of the perturbations 
are comparable with the scale heights of the variations for 
the variables. For k ^ 10~ 10 cm -1 , the propagation speed 
is slightly less than the isothermal sound speed in both 
cases. 

Thus, radiation-modified sound waves can be amplified 
in late-type stars with significant radiative forces. How- 
ever, the perturbations are not very strongly excited with 
minimum growth times ~ 10 5 s, contrary to the findings 
by Carlberg (1980) for hot stars, mainly as a result of the 
much smaller radiative fluxes (F cx T£ s ). For the instabil- 
ity to generate large amplification, the sound waves must 
be reflected at the surface and pass through the region 
with high radiative forces several times. However, the in- 
wards propagating modes are damped for the same reason 
as the outgoing are amplified, with the magnitude of the 
damping rate similar to the growth rate, which suggests 
that the instability will still not be very efficient. In order 
to estimate the behaviour of the instability in the global, 
non-linear regime, a more sophisticated study than pre- 
sented here must be carried out. 



(a) Solar 






s 

/ 

/ 

/ 

/ 

/ 

( 

i 





10 



-10 



10" 



Wavenumber k [cm ] 





(b) RCrB 
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Fig. 5. The numerical solutions to the full set of linearized 
equations for the gravity-acoustic wave modes for a solar and b 
RCrB compositions. The real parts are the dotted curves. The 
solid lines correspond to amplification of the imaginary parts, 
while the dashed curves represent damping. For wavenumbers 
with damping, the absolute magnitude of the imaginary part 
is shown 



5.2. The simplifying case of a isothermal, non-ionizing at- 
mosphere 

In order to gain further theoretical insight into which 

atmospheric properties contribute to the instability pre- De P gas Dp 

sented above, the equations can be simplified by assuming ^ Dt p Dt 



the unperturbed atmosphere to be isothermal and neglect- 
ing ionization and scattering. Furthermore, only the lim- 
its of optically thin and optically thick disturbances will 
be studied and hence the combined radiative energy and 
momentum equation will not be needed. The gas energy 
equation will thus instead be written as: 



= -£ 



(20) 
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where L denotes the net radiative cooling rate of the gas. 
Depending on the optical thickness of the disturbance, C 
will be described by either Newtonian cooling or equilib- 
rium diffusion (e.g. Mihalas & Weibel Mihalas 1984). The 
rate of work done by g r ad has been ignored. 

The linearized continuity equation remains the same 
but the linearized momentum equation reads as: 



"Tn 



ik 



fh 
Po 



where g p = p dg la .d/dp and g T 
energy equation looks like: 



[c\ik + g cB - g p ] + 
-9t] +Vi [-iuj] = 0, (21) 
T dg rad /dT, while the 



fh 
Po 



To 



7- 



1 



C 



(22) 



Here C represents the term arising from the net cooling 
rate C and 7 the ratio of specific heats. For Newtonian 
cooling 



s _ i6/c Ross (T R r^ 

(T-l)cvT 
and in the diffusion regime 
- _ 16K Ross anT$ k 2 



(7 - l)CyT 3k Ross /3 : 



,2 ' 



(23) 



(24) 



C is the inverse of the radiative cooling time-scale 
(e.g. Mihalas & Weibel Mihalas 1984), where the factor 
fc 2 /(3K Ross p 2 ) <C 1 takes into account the optical depth 
effects for optically thick conditions. The simplified dis- 
persion relation then reads: 



-iu 3 + u?l (7 - 1) + 

+iuj [7fc 2 4 + H^ 1 (g cS - g p )] + 

-wk [g p + (7 - 1) g T - g cB + jc^Hp 1 ] + 

-£( 7 -l) [k 2 c 2 T + H; 1 (g cS - g p )] + 

-ikt (7 - 1) [g p - g cS + c^H- 1 ] = 0. (25) 

For k > 10~ 10 cm. - , H^ 1 (g c s — g p ) is much smaller than 
k 2 c\ for the two investigated models and can be ignored. 

In situations of rapid radiative cooling (C 3> w), the 
dispersion relation can be simplified to yield the solutions 



± 



fc CT + — { 3p 



5cff 



^4} 



(26) 



i.e. isothermally propagating, radiation-modified sound 
waves. In the limit of slow cooling (C <C uj), as expected 
adiabatic acoustic waves are recovered: 



± 



kc s + — {g p + (7 - 1).9t - ffeff + H p x c 2 } 



,(27) 



which travel with the adiabatic sound 
speed c s = (dP gas /dp) s = jP gas /p. For k <, 10 _10 cm _1 , 
the propagation speed will be modified by the neglected 
terms H^ 1 (g c ff — g p ), though the growth/damping rates 
will remain the same. 

For both isothermal and adiabatic perturbations, the 
sound waves are amplified by the action of the radiative 
acceleration, as long as the opacity increases with temper- 
ature and density (assuming constant radiative flux): 



9 P 
9T 



c dT ' 



(28) 
(29) 



Furthermore, in a P gas -inversion g cB is negative and hence 
tends to excite an initial disturbance. However, the normal 
momentum conservation counteracts this due to propaga- 
tion in a density stratified medium; the density inversion 
will tend to damp the amplitude for outwards propagat- 
ing waves. In fact, with this simplified atmosphere and the 
values in Table |l|, in both models the last term dominates. 
Hence, the amplitude of sound waves will decay as they 
propagate outwards with these assumptions, despite the 
super-Eddington luminosities and the increase in opacity 
upon com pression. Thus, since the numerical results in 
Sect. 5.1 show amplification for the case without these 
simplifying assumptions, the effects of e.g. ionization and 
a temperature gradient tend to work together with the ra- 
diative force and its derivatives to overcome the damping 
contribution from the density gradient. 

6. Conclusions and speculations 

It has been shown that H-deficient model atmospheres 
may have radiative forces which exceed gravity, simi- 
lar to the case of some luminous H-rich stars (Lamers 
& Fitzpatrick 1988; Gustafsson & Plez 1992). Such 
super-Eddington luminosities manifest themselves as -P gas - 
inversions in hydrostatic atmospheres. The inversions are, 
however, not an artifact of this assumption but can also be 
present in dynamical atmospheres. The inversions are not 
removed unless a very high mass loss rate is present, and 
is neither necessarily removed by convection if the large 
radiative forces occur close to the surface. Thus, super- 
Eddington luminosities does not automatically initiate a 
stellar wind. 

The location of the H-deficient R CrB stars in the im- 
mediate proximity of the computed opacity-modified Ed- 
dington limit (Fig. |^), is certainly suggestive of a connec- 
tion between the enigmatic visual declines of the stars and 
the Eddington limit. It is therefore proposed that insta- 
bilities as the RCrB stars encounter the Eddington limit 
during their evolution towards higher T e g are the unknown 
trigger mechanism for their famous variability. The Ed- 
dington limit may thus be the "smoking gun" for the de- 
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clines: gas is ejected due to the high radiative forces in the 
atmospheres, which then cools rapidly first radiatively and 
later from adiabatic expansion to reach sub-equilibrium 
temperatures and thus possibly enable dust condensation 
(cf. Woitke et al. 1996 for the similar case of shocks instead 
induced by pulsations). Such gas ejections may be observ- 
able as absorption components in strong lines, as observed 
by Rao & Lambert (1997) in RCrB at maximum light. If 
a connection indeed exists, it could explain why the EHc 
and HdC stars do not show R CrB-like variability as a re- 
sult of being located on the stable side of the Eddington 
limit. It also points to an interesting similarity between 
the RCrB stars and the LB Vs. Thus, the two types of 
stars known to be situated at the Eddington limit both 
show eruptive behaviour, which suggests a similar under- 
lying physical explanation for their variability. 

A search for radiative and dynamical instabilities in 
the atmospheres has been carried out for both H-rich and 
H-deficient late-type supergiants, but only partly success- 
fully. Sound waves are found to be amplified by the large 
radiative forces but the linear stability analysis reveal only 
rather slow growth rates compared to for early-type stars 
despite the super-Eddington luminosities. Such radiation- 
modified sound waves may thus be partly responsible 
for the semi-regular pulsational variations of such super- 
giants, but it is doubtful whether the instability is efficient 
enough to eject gas clouds from the atmospheres as specu- 
lated above. The atmospheres are also found to be close to 
dynamically unstable, which might give rise to increased 
mass loss, as previously suggested for LBVs (Stothers & 
Chin 1993) and for the termination of the AGB phase (e.g. 
Paczynski & Ziolkowski 1968; Wagenhuber & Weiss 1994). 
This does not, however, explain in what way the RCrB 
stars are special compared to other similar stars. In fact, 
the H-rich models seem to be more vulnerable to such in- 
stabilities, which suggests that dynamical instabilities are 
not the answer to the declines of the R CrB stars. Another 
possibility is of course that the same dynamical instabili- 
ties occur in both types of stars but that the H-deficient 
and C-rich environment of the R CrB stars greatly favours 
dust condensation, and thus the observational manifesta- 
tion, compared to for H-rich compositions. 

One can also imagine that the Eddington limit is still 
responsible for the behaviour of the RCrB stars though 
in a more indirect way. Langer (1997a, b) has shown that 
the effect of rotation coupled with strong radiative forces 
may lead to increased mass loss in the equatorial plane, 
and it is possible that such a mechanism is at work in 
these supergiants, despite the presumably relatively small 
rotational velocities. The difference in observed variability 
between the R CrB stars and the HdC and EHe stars could 
then be the result of viewing angle. There are some ob- 
servational evidence for such bipolarity (e.g. Rao & Lam- 
bert 1993; Clayton et al. 1997) Violent instabilities due 
to strange mode pulsations leading to dramatically in- 
creased mass loss is another possibility, which is partly 



related to the Eddington limit, since both have their ori- 
gin in ionization zones, as do the dynamical instabilities 
discussed above. Whether any of these instabilities, or a 
combination thereof, could help explain the variability of 
the RCrB stars and the LBVs certainly deserves further 
investigation. 

Finally, it should be remembered that the stability 
analyses presented here assume the unperturbed model 
to be a realistic description of the stellar atmosphere. At 
least for the RCrB stars, there are strong indications that 
this is in fact not the case (Gustafsson & Asplund 1996; 
Asplund et al. 1997a, b,c; Lambert et al. 1997). If the su- 
pergiant atmospheres are indeed distinctly different from 
the predictions of standard models, the efficiency of ra- 
diative instabilities may have been underestimated. One 
can suspect that models based on the normal assumptions, 
such as the mixing length theory for convection, near the 
Eddington limit may be very unsuitable. Clearly, radiative 
hydrodynamical simulations of such atmospheres would 
be valuable, also since they would shed further light on 
possible instabilities. 
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